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ABSTRACT 

We present first results of our planet search program using the 9.2 meter Hobby-Eberly 
Telescope (HET) at McDonald Observatory to detect planets around M-type dwarf stars via 
high-precision radial velocity (RV) measurements. Although more than 100 extrasolar planets 
have been found around solar-type stars of spectral type F to K, there is only a single M-dwarf 
(GJ 876, Delfosse et al. 1998; Marcy et al. 1998; Marcy et al. 2001) known to harbor a planetary 
system. With the current incompleteness of Doppler surveys with respect to M-dwarfs, it is not 
yet possible to decide whether this is due to a fundamental difference in the formation history 
and overall frequency of planetary systems in the low-mass regime of the Hertzsprung-Russell 
diagram, or simply an observational bias. Our HET M-dwarf survey plans to survey 100 M- 
dwarfs in the next 3 to 4 years with the primary goal to answer this question. Here we present 
the results from the first year of the survey which show that our routine RV-precision for M- 
dwarfs is 6 ms~^. We found that GJ 864 and GJ 913 are binary systems with yet undetermined 
periods, while 5 out of 39 M-dwarfs reveal a high RV-scatter and represent candidates for having 
short-periodic planetary companions. For one of them, GJ 436 (rms = 20.6 ms~^), we have 
already obtained follow-up observations but no periodic signal is present in the RV-data. 

Subject headings: stars: late-type — stars: low mass — planetary systems — techniques: radial velocities 



1. Introduction 

M-dwarf stars form the majority of stars in our 
galaxy, both in numbers as well as in total mass. 
However, Doppler surveys looking for stellar re- 
flex motions due to planetary companions via pre- 
cise radial velocity measurements, have tradition- 
ally focused on the brighter F,G and K-type stars 
in order to obtain a sufficient signal-to-noise ra- 
tio of the high resolution (typically R > 50,000) 



^ Based on data collected with the Hobby-Eberly Tele- 
scope, which is operated by McDonald Observatory on 
behalf of The University of Texas at Austin, the Penn- 
sylvania State University, Stanford University, Ludwig- 
Maximilians-Universitat Miinchen, and Georg-August- 
Universitat Gottingen. 



spectra. This has led to the discovery of more 
than 100 extrasolar giant planets orbiting solar- 
type stars (Mayor & Queloz 1995, see e.g. Fischer 
et al. 2003 for recent detections). Due to their 
intrinsic faintness and thus need of large aperture 
telescopes, M dwarfs constitute only a small frac- 
tion of the target-samples of these surveys. Up 
to now, there is only a single M-dwarf, GJ 876 
(M4V, M=0.3 Mq), known to possess a planet 
with a minimum mass of m sin z « 2 Mjup in a 60 
day orbit (Delfosse et al. 1998; Marcy et al. 1998). 
Follow-up observations revealed the presence of a 
second companion with msini — 0.56 Mjup mov- 
ing in a 2:1 mean motion resonance (Marcy et 
al. 2001). In combination with the RV data Bene- 
dict et al. (2002) used the HST Fine Guidance 
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Sensors (FGS) to measure the astrometric pertur- 
bation of GJ 876 caused by the outer planet, and 
thus to determine the orbital inclination and the 
true mass of the companion. 

There were few attempts to target specifically 
M-dwarfs to look for substellar companions. Start- 
ing in 1995 the Delfosse et al. (1999) Doppler sur- 
vey monitors a volume-limited sample of 127 M-- 
dwarfs with the ELODIE spectrograph at Obser- 
vatoire de Haute-Provence and the CORALIE in- 
strument at ESO La Silla. Their long-term RV- 
precision spans from 10 ms~^ for the brighter stars 
to 70 ms"-*^ for the fainter objects. Beside the 
planetary companion to GJ 876, they also dis- 
covered several spectroscopic stellar binary sys- 
tems in their sample (Delfosse et al. 1999). The 
Keck Doppler survey of the Berkeley group in- 
cludes » 150 M-stars and they achieve a long-term 
RV precision of 3-4 ms'^ (Vogt et al. 2000). The 
Keck Hyades program of Cochran et al. (2002) also 
contains a subsample of 20 M-stars. In the south- 
ern hemisphere we have been using the ESO VLT 
and the UVES spectrograph to observe 25 M-stars 
for 2 years now. The high RV precision and sam- 
pling density we obtain for Proxima Cen (M5V) 
would have already allowed us to detect very low- 
mass planets (msini = 0.0126 — 0.019 Mjup or 
msini = 4 — 6 M^) inside the habitable zone of 
that star (Endl et al. 2003), while in the case of 
Barnard's star (M4V) we successfully measured 
the secular acceleration of the RV for the first 
time and discovered a correlation between our RV 
data and the Ha emission level for this M-dwarf 
(Kurster et al. 2003). 

With the exception of the Delfosse et al. and 
our VLT program, none of the above mentioned 
surveys aim for M-dwarfs as primary targets and 
in general low-mass stars constitute only a s\ib- 
sample of larger targetlists. Our HET program 
only observes M-dwarfs, since we are especially in- 
terested in the prevailance of planetary systems in 
this section of the HR-diagram. 

2. Planets around M-dwarfs 

M-dwarfs form the majority of stars in our 
galaxy, both in total numbers as well as in cu- 
mulative mass. Knowledge about the properties 
and frequency of planetary systems for this type 
of stars would therefore have vast implications on 



the overall statistics of planetary companions in 
the galaxy. 

So far no detailed physical models for planet 
formation in general or as a function of stellar 
mass in particular are available. To our knowl- 
edge, Nakano (1988) was the first to study planet 
formation timescales for stars of varying mass and 
lifetimes. Wetherill (1996) simulated the forma- 
tion of terrestrial planets for stars of different mass 
and disks of different surface densities, and found 
that the end-products of planetesimal accretion 
are rather insensitive to the stellar mass. How- 
ever, the low-mass star chosen in both studies has 
a mass of 0.5 M^, which is the high end of the 
M-dwarf mass distribution. Since formation of gas 
giant planets is even less understood (it is not clear 
yet whether jovian planets form in a "bottom-up" 
process, i.e. by gas accretion onto a rocky core 
of w 10 Me (e.g. Pollack et al. 1996), or by a 
type of gravitational disk instabilities leading to 
a "top-down" collapse and rapid formation of the 
gas giant without the need of a rocky core (e.g. 
Boss 1997), reliable theoratical predictions on the 
frequency of planetary systems depending on ini- 
tial conditions like stellar primary and disk mass 
as well as disk composition are difficult to make. 

Therefore, our goal is to obtain a useful statisti- 
cal overview of planetary systems in the M-dwarf 
regime by observational means. In the near future 
we will know if GJ 876 remains a special case or 
if more planets around M-dwarfs will be detected 
by our survey or one of the other programs. 

3. The HET survey and radial velocity re- 
sults 

Table 1 lists 35 targets (out of a total of 39, 
4 candidate objects for short-periodic companions 
are not listed) of our survey, the total rms-scatter 
of their RV-data and the duration of monitoring. 
With a few exceptions the short time span of ob- 
servations limits our detection capability to short- 
periodic planets with periods of a few days. Ex- 
tension of the time baseline will allow us to become 
sensitive also to longer period planets at larger or- 
bital separations. 

All targets were selected based on the Gliese 
catalogue of nearby stars (Gliese & Jahreifi 1991) 
and the Hipparcos astrometric database (ESA 1997) 
In order to minimize the effects of stellar activity 
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on our RV measurements (especially rotational 
modulation by star spots) we choose only tar- 
gets with no detected or low X-ray emission using 
the ROSAT all-sky survey results from Hiinsch et 
al. (1999). 

In our observing strategy we usually first take 
5 RV measurements per target over a short time 
to test for short-periodic variability. If variablity 
is detected the same target will be re-scheduled 
for more intense monitoring. Seemingly constant 
stars will be re-observed later in our survey and 
with a longer cadence in order to check for long- 
term variability (companions at larger orbital sep- 
arations). 

This work presents the first results our pro- 
gram, which we plan to continue until we have 
sampled « 100 M-dwarfs with the necessary RV- 
precision to find orbiting giant and lower mass gi- 
ant planets. 

For the observations wc use the HET High- 
Resolution- Spectrograph (HRS; Tull 1998) which 
is located in an insulated chamber below the tele- 
scope. An optical fiber connects the instrument to 
the prime focus of the HET and also provides ex- 
cellent scrambling of the light beam in order to iso- 
late the spectrograph from variations in the tele- 
scope pupil illumination. All spectra of our pro- 
gram are taken at a resolving power of i? = 60, 000 
and an I2 cell is inserted into the light path for 
self-calibration purposes. 

3.1. RV-precision 

To extract the Doppler shift information from 
the M-dwarf spectra with the superimposed io- 
dine vapor (I2) reference spectrum we perform a 
full data modeling including the reconstruction of 
asymmetries of the instrumental profile (IP) at the 
time of observation. In this modeling scheme we 
subdivide the spectrum into smaller segments and 
establish a best-fit model, based on a stellar tem- 
plate spectrum of the target and a laboratory spec- 
trum (a Fourier Transform Spectrometer scan) of 
the I2 cell, for each segment. The RV of the target 
is measured differentially to the template while the 
dense forest of superimposed I2 lines delivers the 
wavelength calibration and the information to re- 
construct the IP. The uncertainty of the RV mea- 
surement is determined from the overall scatter 
of all segments along a spectrum. A detailed de- 



scription of this technique can be found in Butler 
et al. (1996) and Endl, Klirster, Els (2000). 

Fig. 1 shows the histogram of our RV-results. 
The bulk of the distribution on the left side con- 
sists of the M-dwarfs, which are either constant 
within our measurement precision or their resid- 
uals after subtraction of RV-trends do not reveal 
additional scatter (GJ 623, GJ 708, GJ 864 and 
GJ 913). Fitting a Gaussian to this part of the 
histogram yields a mean value of 5.7 ms"-"^ with 
a = 2.2 ms~^, which is our present mid-term 
RV-precision for stars in the magnitude range of 
F = 9 - 11 mag (with the bright {V = 7.5) M2V 
star GJ 411 being the sole exception). The right- 
hand "tail" of the RV-scatter distribution con- 
tains low-amplitude variable stars (all with rms 
> 10 ms~^), which we want to sample more fre- 
quently to check for short-periodic planetary com- 
panions as possible cause of the higher RV-scatter. 

There are several sources which contribute to 
the RV-distribution seen in Fig. 1: 1. intrinsic 
variability of the targets (i.e. activity induced 
"RV-jitter" or Keplerian motion due to unknown 
substellar or stellar companions), 2. the photon 
noise, 3. the stability of the spectrograph used 
for the RV measurements, and 4. "algorithmic 
noise" introduced during the data modeling pro- 
cess (e.g. inadequate characterization of the in- 
strumental profile and imperfect stellar template 
deconvolution) . For the M-dwarf HET sample we 
conclude that photon noise is the major contrib- 
utor to the RV-scatter. For earlier type stars of 
visual magnitude of F w 7 mag, using the same 
intrumental setup and analysis technique, we ob- 
tain a routine long-term RV-precision of « 3 m s""'^ 
(Cochran et al. 2003), a factor 2 better than for 
the fainter M-dwarfs. 

3.2. RV-results 

The RV-results of all targets listed in table 1 are 
displayed in Fig. 2 to Fig. 10. To allow a better 
comparison we plot all results on the same velocity 
scale (—30 ms~^ to -1-30 ms~^), with the excep- 
tion of GJ 436 (Fig. 6) where the scale is enlarged. 
For GJ 623, GJ 748, GJ 864 and GJ 913 we show 
the RV-data after subtraction of a large amplitude 
RV trend (see section 3.3). 

For one of the targets in the high-rms part of 
the distribution, GJ 436, we have already obtained 
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follow-up observations (sec Fig. 6). A search 
for periodic signals using the Lomb-Scargle peri- 
odogram (Lomb 1976, Scargle 1982) did not reveal 
any significant signal in the power spectrum (the 
highest peak has a false-alarm-probability of 32%). 
The X-ray luminosity of Lx = 0.7 x 10^'^ ergs~^ 
(Hiinsch et al. 1999) might point toward stel- 
lar activity as cause of the excess scatter. How- 
ever, GJ 411 with a similar level of coronal emis- 
sion (Lx = 0.6 X lO^"^ crgs^i) and GJ 272 with 
Lx = 3.1 X 10^'' ergs~^ the brightest X-ray source 
in our sample, do not show an increased RV- 
scattcr. The cause of the variability of GJ 436 
hence remains unknown. The other variable tar- 
gets are currently in the observing queue of the 
HET and are pending re-observations. 

3.3. Linear and non-linear RV trends 

Three targets of our sample are either pre- 
viously known (GJ 623, GJ 748) or suspected 
(GJ 913) binary stars. 

GJ 623 was found to be an astrometric binary 
by Lippincott & Borgman (1978) with the latest 
published value for the orbital period of 3.7 yrs 
(Nidever et al. 2002). Our HET RV-results are 
shown in Fig. 11 along with a parabolic fit to 
the data. Fig. 7 displays the RV-residuals after 
subtraction of this trend. This data-set is now 
being used to further improve the binary orbital 
solution and, in combination with HST FGS as- 
trometry, to accurately determine the masses of 
the components of the GJ 623 system (Benedict 
priv. communication) . 

G J 748 is another astrometric binary star (Har- 
rington 1977), where we only see a small fraction 
of the P = 2.45 yrs (Franz et al. 1998) orbit 
(Fig. 12). The residual scatter around a parabolic 
fit is 5.5 ms"^ (Fig. 9). 

The entry of GJ 913 in the Hipparcos catalogue 
contains a double/multiple system annex flag X 
(stochastic model), which is a hint that the star 
is actually a short-periodic binary. The HET RV- 
data (Fig. 13) reveal a large amplitude variation 
with a Dopplershift of ~ 500 ms~^ over the time 
span of 20 days, indicative of binary orbital motion 
with an unknown period and thus confirming the 
Hipparcos results. Further observations will allow 
us to determine a spectroscopic orbit for this star. 

In the case of GJ 864 we discovered a linear ac- 



celeration in the first 5 RV measurements. Since 
GJ 864 is not a known (or suspected) binary star 
and was classified as an RV-constant star in the 
CORAVEL survey (Tokovinin 1992), we immedi- 
ately re-scheduled GJ 864 for follow-up observa- 
tions. As can be seen in Fig. 14 the linear RV- 
trcnd continued until the star became unobserv- 
able by the HET in December 2002. Using aux- 
iliary RV measurements with the McDonald 2.7 
meter telescope in January 2003 we found that 
the RV-trend of GJ 864 remained linear without 
any sign of curvature and thus the best explana- 
tion for the RV-data is a stellar companion at yet 
undetermined separation. 

4. Discussion 

So far, we have found 5 stars out of a sample of 
39 M-dwarfs which show a higher RV-scatter than 
the rest of the sample and which we plan to mon- 
itor intensively in the near future to check if the 
cause of the variability are short-periodic plane- 
tary companions. For GJ 436, the one target for 
which we already obtained follow-up observations, 
we could not find a periodic signal in the RV-data 
and the cause of variability remains unknown. 

Two stars of our sample turned out to be previ- 
ously unknown binary stars: GJ 864 and GJ 913. 
For both targets we soon might be able to estab- 
lish a spectroscopic orbital solution. 

Despite their faintness M-dwarfs hold the 
promise to detect lower mass extrasolar planets 
than what is usually possible by precise Doppler 
surveys. Due to the lower stellar primary mass 
even planets with the mass of Neptune or lower 
become detectable in short-periodic orbits. For 
instance, a planet with msini = 17 Mq in a 4 day 
orbit induces an RV semi-amplitude of ~ 11 ms~^ 
on an MOV (0.5 Mq) star, a 2f7-signal for our 
program. For later M-dwarfs these detection 
limits approach the terrestrial planet regime as 
demonstrated by the mass-limits for Proxima Cen 
(Endl et al. 2003) and Barnard's star (Kiirster 
et al. 2003). Thus, for the time being only pre- 
cise Doppler monitoring of these low-mass stars 
has the capability to detect planets with masses 
below 10 M0. Supposendly, such low-mass plan- 
ets would have formed by planetesimal accretion, 
very similar to the way the Earth was formed, 
rather than by the processes which are responsi- 
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ble for the formation of gas giant planets (of a 
few hundred M®). Clearly the detection of these 
kind of planets would be of special interest for 
future space missions which will be dedicated to 
search for Earth analogs (KEPLER, NASA's TPF 
& ESA's Darwin). 

After the first year of the project the size of 
our target sample is still insufficient to arrive at 
any significant statistical conclusions, but we con- 
stantly add new targets to onr list and the exten- 
sion of our survey, both in target numbers as well 
as in time, will allow us in the future to deter- 
mine useful quantitative limits for planets around 
M-dwarf stars. 
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Table 1: Radial velocity results for 35 M-dwarfs (4 
variable M-dwarfs of our program are not listed 
here). Spectral classification and visual magni- 
tudes are taken from the Gliese catalogue, N is 
the total number of obtained RV measurements, 
rms is the total scatter of the RVs and AT de- 
notes the duration of monitoring (i.e. the time 
span from first to last observation of this star). 
Notes: for the 4 stars marked with an asteriks the 
rms is given after subtraction of an RV-trend (see 
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Fig. 1.— Histogram of the HET/HRS RV-results 
for 39 M-dwarfs from our sample. The left-hand 
side contains constant stars and targets which 
show no excess residual signal after subtraction of 
RV-trends. Fitting a gaussian to this part of the 
histogram yields a mean RV-scattcr of 5.7 ms~^ 
and a = 2.2 ms^^. The right-hand "tail" of the 
distribution consists of more variable stars with 
rms- values higher than 10 ms~^. 
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Fig. 2.— RV results for GJ 2, GJ 87, GJ 155.1 
and GJ 184. 
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Fig. 3.— RV results for GJ 192, GJ 
GJ 251.1, and GJ 272. 



3352, Fig. 4.— RV results for GJ 277.1, GJ 281, GJ 289, 

and GJ 308.1. 



8 



GJ 310 



RMS = 15.9 m/s 



53660 52680 52700 

JD 2.400,000+ [days] 



GJ 2085 



RMS = 6.2 m/s 



52660 52662 52664 52666 52668 

JD 2.400,000+ [days] 



GJ 328 



RMS = 6.2 m/s 



52640 52650 52660 52670 

JD 2.400.000+ [days] 



GJ 436 




RMS = 20.6 m/s 



52300 52400 52500 52600 

JD 2.400.000+ [days] 



GJ 378 



RMS = 4.6 m/s 



52650 5S652 53654 

JD 2.400.000+ [days] 



GJ 535 



RMS = 3.0 m/s 



52400 5^4 10 

JD 3,400.000+ [days] 



GJ 411 



RMS = 7.2 m/s 



52300 52400 53500 

JD 3.400,000+ [days] 



GJ 552 



RMS = 11.3 m/s 

52450 



52460 52470 
JD 3.400,000+ [days] 



Fig. 5.— RV results for GJ 310, GJ 328, GJ 353, 
GJ 378, and GJ 411 



Fig. 6.— RV results for GJ 208-5, GJ 436 (note 
the larger velocity scale), GJ 535, and GJ 552. 



9 



GJ 563.1 



RMS = 4.9 m/s 



S2420 52425 52430 52435 52440 52445 

JD 2.400.000+ [days] 



GJ 699 



RMS = 7.2 m/s 



52420 52430 52440 52450 

JD 2.400.000+ [days] 



GJ 623 residuals 



RMS = 5.0 m/s 



53500 
JD 3,400,000+ [days] 



GJ 708 



RMS = 9.2 m/s 



53435 53430 53435 53440 53445 

JD 3,400,000+ [days] 



GJ 671 



RMS = 6.1 m/s 



52446 52448 52450 52453 52454 52456 

JD 2,400.000+ [days] 



GJ 4070 



RMS = 8.0 m/s 

52400 



52430 52440 
JD 3,400.000+ [days] 



GJ 687 






i - 


RMS = 2.9 m/s 





53398 52400 53403 

JD 3.400,000+ [days] 



GJ 731 



{ + 



RMS = 3.4 m/s 



53444 53446 53448 

JD 3.400,000+ [days] 



Fig. 7.— RV results for GJ 563.1, GJ 623 (resid- 
uals), GJ 671, and GJ 687 



Fig. 8.— RV results for GJ 699, GJ 708, GJ 4070, 
and GJ 731. 
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Fig. 9.— RV results for G.J 748 (residuals), 
GJ 839, GJ 864 (residuals), and GJ 849. 
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Fig. 10.— RV results for GJ 895 and GJ 913 
(residuals). 
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Fig. 11.— HET RV-data for the known binary 

GJ 623 showing a part of the binary orbit. The 
rms scatter around the best-fit parabolic trend is 
5.0 ms~^ (see Fig. 7 for the residuals). 
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Fig. 12.— HET RV-data for the binary GJ 748 

showing a small fraction of the binary orbit [P = 
2.45 yrs). The rms scatter around the best-fit 
parabolic trend is 5.5 ms~^ (see Fig. 8 for the 
residuals). 
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Fig. 14.— HET RV-data of GJ 864 showing a 
linear acceleration of —5.5 ms^^d^^. The rms- 
scatter around this trend is 7.3 ms^^. Using addi- 
tional observations with the McDonald 2.7 meter 
telescope we confirmed that the linearity of this 
trend continued until January 2003 and that the 
best explanation is a previously unknown stellar 
companion at yet undetermined separation. 
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Fig. 13. — GJ 913, a Hipparcos candidate for a 
short-period binary system, shows a large ampli- 
tude RV variation indicative of a stellar compan- 
ion. The rms scatter around the best-fit parabolic 
trend is 6.5 ms~^ (see Fig. 9 for the residuals). 
Further observations might allow us to determine 
a spectroscopic orbital solution. 
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